Astronomy & Astrophysics manuscript no. 10969AA 


© ESO 2008 


November 6, 2008 





Studying the properties of the radio emitter in LS 5039 

(Research Note) 



V. Bosch-Ramon 



OO 

O 

o 

> 
O 



O 



> 
00 

O 



OO 
O 



Max Planck Institut fiir Kemphysik, Saupfercheckweg 1, Heidelberg 691 17, Germany; vbosch@mpi-hd.mpg.de 
Preprint online version: November 6, 2008 

ABSTRACT 

Context. LS 5039 is an X-ray binary that presents non-thermal radio emission. The radiation at ~ 5 GHz is quite steady and optically 
thin, consisting on a dominant core plus an extended jet-like structure. There is a spectral turnover around 1 GHz, and evidence of 
variability at timescales of 1 yr at 234 MHz. 

Aims. We investigate the radio emitter properties using the available broadband radio data, and assuming two possible scenarios to 
explain the turnover: free-free absorption in the stellar wind, or synchrotron self-absorption. 

Methods. We use the relationships between the turnover frequency, the stellar wind density, the emitter location, size and magnetic 
field, and the Lorentz factor of the emitting electrons, as well as a reasonable assumption on the energy budget, to infer the properties 
of the low-frequency radio emitter Also, we put this information in context with the broadband radio data. 

Results. The location and size of the low-frequency radio emitter can be restricted to ^ few AU from the primary star, its magnetic 
field to ~ 3 X 10"^ - 1 G, and the electron Lorentz factors to ~ 10 - 100. The observed variability of the extended structures seen 
with VLB A would point to electron bulk velocities ^ 3 x 10^ cm s"', whereas much less variable radiation at 5 GHz would indicate 
velocities for the VLBA core S 10** cm s"'. The emission at 234 MHz in the high state would mostly come from a region larger than 
the dominant broadband radio emitter 

Conclusions. We suggest a scenario in which secondary pairs, created via gamma-ray absorption and moving in the stellar wind, are 
behind the steady broadband radio core, whereas the resolved jet-like radio emission would come from a coUimated, faster, outflow. 
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1. Introduction 

■ LS 5039 is an X-ray binary that presents broadband non-thermal 
emission: in radio (e.g. Marti et aL 1 19981 M98 hereafter); in X- 
rays (e.g. Bosch-Ramon et al. l2007l l; at GeV energies (Paredes et 

■ al.gOOOi PAO); and at TeV energies (Aharonian 2005, 2006), as 
I well as extended ~ 10 - 1000 AU radio jets (e.g. PAO; Paredes 
. et aL 120021 Ribo et al. |2008l -R08 hereafter-). The mass of the 
' compact object is not well constrained (Casares et al. 120051 

C05 hereafter), and it may harbor an accreting black-hole (e.g. 
Paredes, Bosch-Ramon & Romero |2006l ) or a non-accreting pul- 
sar (e.g. Martocchia, Motch & Negueruela|2005i Dubus |2006t . 
The source is at a distance of 2.5 kpc (C05). 

Low-frequency (^ 1 GHz) radio emission has been detected 
by the Giant Metrewave Radio Telescope -GMRT- (Pandey et al. 
120071 Godambe et al. 120051 P07 and G08 hereafter) in LS 5039. 
This radio emission is very likely of synchrotron origin, match- 
ing well with the non-thermal radio spectrum at higher fre- 
quencies (see fig. 2 in G08). G08 found a spectral turnover at 
Vmax = 964 ± 104 MHz (with a flux F^^^ ^ 40 mJy), fluxes 
of ^ 17.0 + 1.1 and 34.2 + 2.8 mJy at 234 and 614 MHz, re- 
spectively, and a spectral index a = 0.75 + 0.11 (Fy oc v") 
in this frequency range. These data were taken simultaneously 
(MJD=53788.78). Above the spectral turnover, at 1.28 GHz, the 
flux was 38.1 + 2.3 mJy, being this data point obtained 18 days 
later (MJD=53806.73) than those at lower frequencies. Above 
the spectral turnover, a was » -0.43, very similar to that found 



by M98. The P07 observations, taken in two different epochs 
(MJD=53220 and 53377) more than 1 yr apart from the G08 
observations, found significantly different fluxes at 234 MHz: 
65. 04 ±7. 71 and 74.56 + 5.63 mJy, which is evidence of variabil- 
ity. Otherwise, the fluxes of P07 at 614 MHz were 28.55 + 3.56 
and 34.29 + 2.34 mJy (MJD=53220, 53377), similar to the value 
of G08. The a found by P07 were ^ -0.81 and -0.86. 

The spectrum found at low frequencies by G08, far more 
inverted than the one coiTesponding to a monoenergetic distri- 
bution of particles (i.e. a - 1/3), could not be explained by 
an electron population with a high minimum energy, being the 
most likely origin of the turnover either free-free absorption in 
the stellar wind or synchrotron self-absorption. Based on this, we 
investigate in this work the properties of the low-frequency ra- 
dio emitter in LS 5039 adopting a plain framework in which the 
radio emission, at different naiTow frequency bands, is produced 
in an homogeneous spherical emitter with magnetic field B, size 
r, and at a distance R ^ r from the star/binary system. We also 
derive physical constraints of the broadband radio emitter from 
data at higher frequencies. Finally, we suggest a scenario for the 
radio emission, and a more complex model is applied to test its 
feasibility. 



2. Studying tlie low-frequency radio emission from 
LS 5039 
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In this section, we infer basic properties of the low-frequency ra- 
dio emitter in LS 5039 accounting for the detection of a spectral 
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turnover at w 1 GHz and the observational evidence of variabil- 
ity at 234 MHz. 

For illustrative purposes, we present in Figure [T] a sketch of 
the scenario considered here. In the picture, the low-frequency 
radio emitter location is illustrated for two cases: a low- 
frequency turnover produced by free-free absorption, with the 
emission coming from a region at few AU from the star, or by 
synchrotron self-absorption, with the emitter located at simi- 
lar or larger distances (see below). The wind region optically 
thick by free-free absorption to 1 GHz emission is shown. The 
core and the jet-like emission at 5 GHz, discussed later, are also 
sketched. 

2. 1 . Free-free absorption 

The strongly inverted low-frequency radio spectrum could be 
produced by free-free absorption in the stellar wind. The dis- 
tance at which the stellar wind becomes optically thin to ra- 
dio emission can be computed fixing the free-free opacity Ty fj 
(Rybicki & Lightman [T979] l to 1: 

T. ff ^ 30M^_,v-^^^ ^ax^Au^wlaT'wT = 1 ' (1) 
which, for vghz max = Vmax/1 GHz = 1, yields: 
R«3M%'v-JllT-J{^A\i, (2) 

where M-q - yr"' is the stellar mass loss rate, 

/^Au = Rll AU, Vw 8.3 = Vw/2 X 10'^ cm s"' the stellar wind 
velocity, and 4 = Tw/lO^ K the wind temperature. Adopting 
Jw 4 = 1 (Ki-tidka & Kubat'SOOT), V^j = 1 (COS), and M-7 = 1, 
we obtain » 3 AU (4.5 x 10'^ cm; recall that r $ /?) as the 
distance between the low-frequency radio emitter and the star. 

A M somewhat lower than the values given by COS has 
been taken to account for possible wind dumpiness in LS S039 
(Bosch-Ramon et al. l2007l l. Actually, wind dumpiness may pre- 
vent the wind to be an effective free-free absorber for large 
enough porosity lengths (e.g. Owocki & Cohen |20061 l. In such a 
case, if a pulsar would be present in the system, the radio pulse 
may not be smeared out, hence being detectable. No radio pul- 
sations have been found up to now. 

If the spectral turnover is produced by wind free-free ab- 
sorption, the synchrotron self-absorption frequency must be < 
1 GHz. Thus, to find a restriction for B, one can use the follow- 
ing relationship between Vmax, ^max, and B in the synchrotron 
self-absorption case: 

V a; 3 9 X lO'^F^^^ o 1/5 -4/5 QTT /o^, 
Vmax~-5-VXiU max mJy'^G '^cm ■ ) 

where f ^ax mjy = ^'max/mJy, and r^m = r/cm. Since f max mjy = 
^max/mJy = 40 and r % R ^ 4.S x lO'-' cm, the condition Vmax < 
1 GHz imposes that B < 3 x 10"^ G. 

The constraint for B, together with the following formula 
(e.g. Pacholczvk il970i : 

v«4.2x lO^^Bcr^GHz, (4) 

allows us to constrain the Lorentz factor of the electrons emitting 
at 1 GHz to y > lOl 

Further restrictions on r, B and y can be derived from en- 
ergy budget constraints. On one hand, the synchrotron cooling 
timescale for 1 GHz emitting electrons is: 

fsync ~ S X l(fB-^y-^ ^ 3.3 X lO^B^ s . (S) 

In general, radio electrons will be transported/advected out of 
the emitter well before radiating all their energy via synchrotron 



emission, i.e. fadv ^ fsync- For reasonable advection velocities 
Vadv S 2 X 10^ cm s"' and emitter sizes r S ^ 4.S x lO'^ cm 
(free-free absorption case), one gets fadv ^ 2 x 10^ s, where 
v'adv has been taken assuming that the emitter flow cannot move 
slower than the stellar wind. On the other hand, the total broad- 
band non-thermal luminosity of the source, ~ lO"'* erg s"', 
can be taken as an upper-limit for the injected luminosity in 
the form of radio emitting electrons. Finally, the observed low- 
frequency luminosity is Lghz ~ lO-'" erg s"'. From all this, we 
can write: 

fsync ^ fadv (ie/^GHz) ^ 2 X lO'^ S . (6) 

Therefore, accounting for Eqs. (|4]i, Q and (|6]l, we derive B 
3 X 10"^ G and y ^ 3 x 10^, which allow us to restrict B to 
~ (3 - 30) X 10-^ G, and y to ~ (1 - 3) x 10^. The minimum 
B-value, and v^ax < 1 GHz, imply that r ~ (3 - 4.S) x 10'^^ cm. 

2.2. Synchrotron self-absorption 

If synchrotron self-absorption is the origin of the turnover, we 
can obtain, from Eqs. ([3]) and (|4]i, v^ax = 1 GHz, and Fmax ~ 
40 mJy, the following r - B and r - y relationships: 

r*10'X^'*4xlO'V'''cm. (7) 

Since the synchrotron emitting electrons must be relativistic, y ^ 
10, and therefore, r ^ 10'"^ cm and B < 1 G. 

In this case, having a constraint on r, the argument of the 
limited energy budget can be applied as well, allowing us to 
derive further restrictions on the parameters: B ~ 10"^ - 1 G; 
y ~ 10 - 10^; and r ~ (4.S - 10) x 10'^ cm ~ few AU (recall 
R > r). 

2.3. Variable emission at 234 MHz 

The combination of P07 and G08 observations show variabil- 
ity at 234 MHz of a 400% at year timescales, too large to be 
attributed only to interstellar medium scintillation effects. The 
flux at 614 MHz does not show significant changes. The value 
of fsync at 234 MHz is: 

fsync = 1.4 xlO'^Bc'^'s, (8) 

probably too long to explain the variability just with changes of 
the electron injection. Some level of variability around 1 GHz 
may be induced by changes in the stellar wind if free-free ab- 
sorption produced the spectral turnover, but at the estimated 
distance for the free-free absorption case, synchrotron self- 
absorption should still affect the radiation at 234 MHz due to 
the region compactness. 

The optically thin nature of P07 data, and the fast variability, 
can be otherwise explained by advection with timescales shorter 
than ~ 1 yr together with injection variations. Taking Vadv ~ 
10'" cm s"', ~ 1 yr timescale restricts the size where the (high- 
state) 234 MHz radiation is produced to r 5 3 x 10'^ cm. Vadv 
has been taken here to be at most mildly relativistic to avoid 
the complexities of significant boosting effects (otherwise not 
expected; see Paredes et al. 120001 |2002| l. In addition, the fact 
that this emission is not synchrotron self-absorbed in the P07 
data, and a flux at 234 MHz of 70 mJy, implies a size: 

r>9xlO"*Bc''^~4xlO'V'^^cm. (9) 

Accounting for energy budget constraints, the 1 yr variabil- 
ity, and the constraint y ^ 10 (i.e. electrons must be relativistic). 
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Fig. 1. Sketch of the regions producing radio emission in LS 5039. The main components of the radio emitter at low and high 
frequencies are shown, together with the relevant spatial scales. 



to Eq. we get: B ~ 4x 10* - 0.6 G; y ~ 10 - 400; and 
r ~ R ~ 10'"* - 3 X lO'^ cm. R has similar limits to those of r; 
if bigger, the source would have appeared as extended under the 
GMRT resolution of ~ 10 arc-seconds at 234 MHz. 

Concluding, the high-state 234 GHz emitting region proba- 
bly has a larger size and smaller magnetic field than the region 
producing the emission > 614 MHz and at 234 MHz in the low- 
state, which indicates that they are probably different. The differ- 
ent spectral shape of P07 data and data > 614 MHz would also 
point to a different electron population. Otherwise, despite the 
variability at 234 MHz, the quite steady 614 MHz flux makes 
the conclusions of Sect. 12.1! still valid concerning the break at 
^ 1 GHz. 

In Table [T] at the top, the results of all this section are sum- 
marized. 



3. The broadband radio emission 

3.1. The 5 GHz VLBA radio emission 

The core seen with VLBA at 5 GHz in LS 5039 is not point 
like, presenting a projected size of ~ 8 AU, and two additional 
south-east (SE) and north-west (NW) components at a projected 
distance of ~ 8 AU from the core center are also detected (PAO; 
R08). The variation of the core flux at 5 GHz within 5 days of 
difference is about a 10%, whereas the changes in flux of the 
SE and NW components are of about a 44 and a 260%, respec- 
tively (R08). Assuming that the flux changes of the SE and NW 
component fluxes are showing power engine variations (unnec- 
essarily affecting both components in the same way), plus the 



fact that the synchrotron cooling timescales of radio electrons 
are too long to explain such a variability, we can derive that 
Vadv ^10** cm s"' to have a steady core, and ^ 3 x 10*^ cm s"' 
for the variable SE and NW components. If shorter time vari- 
ability in the SE and NW components were observed, it would 
increase the Vajv lower limit to 3 x 10*'(f/5 days)"' cm s"'. The 
main properties of the VLBA radio emission are summarized in 
Tabled bottom. 

Since LS 5039 harbors a very bright star (C05) and produces 
TeV emission (Aharonian et al. 2005), secondary pairs can be ef- 
ficiently produced in the stellar wind via gamma-ray absorption. 
This population of secondary pairs could be behind the radio 
core (e.g. Bosch-Ramon, Khangulyan & Aharonian l2008at BOS 
hereafter), moving collectively with the stellar wind, hence the 
low inferred advection speed. The feasibility of this scenario is 
illustrated in Fig.|2] where we show the computed radio spectral 
energy distribution (SED) of the secondary synchrotron emis- 
sion created by gamma-ray absorption in LS 5039. The mag- 
netic field in the stellar surface has been fixed to 250 G (expected 
to be reasonable; see Bosch-Ramon, Khangulyan & Aharonian 
2008b), the injected TeV luminosity to 4 x 10^^ erg s"' (simi- 
lar to that inferred by Khangulyan, Aharonian & Bosch-Ramon 
2008 ), and the TeV emitter distance to the stellar companion to 
3.5 X 10'^ cm, not exactly inside the binary system (as proposed 
in Bosch-Ramon et al. 2008b) . 

As seen in the Fig. |2l the computed fluxes match the ob- 
served values, and most of the 5 GHz emission is predicted to 
come from a region of few AU size, similar to first order to 
the values discussed in the previous paragraph. The computed a 
around 5 GHz is ^ -0.4, roughly the same as the value found by 
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Table 1. Basic physical parameters of the radio emitter in LS 5039 
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Fig. 2. Computed radio SEDs of the secondary synchrotron 
emission created by gamma-ray absorption in LS 5039. The 
SEDs corresponding to the total emission, and to that produced 
at distances > 1 AU, from the star are shown. The total differen- 
tial flux at 5 GHz would be w 20 mJy and, > 1 AU, ~ 10 mJy. 
The observed luminosities from the 5 GHz core, and at 234 MHz 
(low state), 614 MHz and 1.28 GHz, are also shown (PAO; R08; 
GOB). 



M98 for the VLA radio emission (which should be dominated by 
the VLBA core). The computed low-frequency spectrum, similar 
to the observed one (G08), suggests that this radiation may come 
also from secondary pairs, although we note that synchrotron 
self-absorption has been roughly computed supposing an homo- 
geneous emitter as wide as long. The magnetic field where most 
of the emission < 1 GHz is produced, at a distance of several AU 
from the star, is ~ 1 G, in agreement with the values inferred in 
Sect. |2] and also compatible with those in the wind at few AU 
from a massive star (Benaglia & Romero '2003'). 

Concerning the SE and NW components, the larger Vadv 
would point to a diff'erent physical structure than that producing 
the radio core, perhaps a real jet. Small changes in the orientation 
of » 12° of this jet-like structure have been found (ROB). This 
might be explained, for instance, by an asymmetry introduced 
by the secondary pairs transported in the stellar wind (BOS), the 
propagation of which might not be spherical, or by jet-stellar 
wind interactions (Perucho & Bosch-Ramon |200^ . 



3.2. The jet-like structure and the 234 MHz variability 

About a 30% of the radiation at 5 GHz comes from a jet-like 
structure of projected size ~ 10 - 1000 AU detected with VLBA, 
MERLIN and EVN (PAO; Paredes et al. 12002! ROB). These ex- 
tended structures may be associated with the 234 MHz high- 
state emission, due to the size constraints for this region inferred 
above. 

The variability seen at 234 MHz may be related to changes in 
the outflow properties (e.g. via hydrodynamical instabilities), but 
not to long-term variations in the primary power supply, since 
it should aff'ect as well the core radio emission. In this frame- 
work, the 234 MHz high state should have associated flux varia- 
tions at higher frequencies visible outside the radio core. Given 
that fadv depends on the size and Vadv, the MERLIN and EVN 
~ 100 - 1000 AU structures at 5 GHz would show variability at 
timescales of at least months. This 234 MHz/large scale jet link 
can be tested with simultaneous radio VLBI observations at low 
and high frequencies taken in two different epochs. We note that 
a radio emitter of r ~ 1000 AU varying with ~ 1 yr timescales 
would have associated a v'adv ^ 10^ cm s-', compatible with the 
Vadv derived for the VLBA SE and NW components. 

4. Discussion 

A plausible scenario for the radio emission in LS 5039 would be 
radiation from secondary pairs created via gamma-ray absorp- 
tion (see BOS), for the compact emitter or core {R ~ few AU), 
plus a coUimated jet, for the resolved emission. 

The jet radio emission from the smaller spatial scales would 
be probably flat, as typically expected in low-hard state micro- 
quasars (e.g. Fender 170011) . but the optically thin radiation from 
secondary pairs would overcome this component, at least up to 
few AU. This would explain why the radio spectrum of LS 5039 
does not look flat (M9S). As observed, this compact radio emis- 
sion would not present significant changes along the orbit, pro- 
vided first that the system is not particularly eccentric and, sec- 
ond, that secondary pairs emit radio emission in a region big 
enough (~ few AU) to smooth out any possible modulation of 
the primary gamma-ray injection. On the other hand, the VLBA, 
MERLIN and EVN extended radio emission would come from 
collimated and steady jets. In short, most of the ^614 MHz 
and the low-state 234 MHz emission would come from the sec- 
ondary pair radiation, forming the radio core, whereas the high- 
state 234 MHz emission, and the 30% of the whole 5 GHz radi- 
ation, would come from the extended jets. 

The 10-1000 AU collimated jets would be consistent with a 
microquasar nature for LS 5039, and hard to reconcile with the 
pulsar cometary-tail model (Dubus |2006t Dhawan et al. I2006I I. 
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although further theoretical work is required before discarding 
out the pulsar scenario, given the complexities inherent to hydro- 
dynamical and magnetohydrodynamical flows (for recent simu- 
lations on pulsar winds, see, e.g., Romero et al. |2007t Bogovalov 
et al. 12008). 

The scenario proposed here can be tested looking for: vari- 
ability in the extended emission at scales ~ 10 AU (during which 
the radio core should remain quite steady); further evidence of 
persistent radio jets beyond the core at scales ~ 10 - 1000 AU; 
hints of a spiral-like radio structure (BOS) in the still marginally 
resolved radio core; some kind of coiTelation between the jet- 
like radio emission at distances < 100 AU and the TeV radia- 
tion (which varies periodically along the orbit; Aharonian et al. 
I2006ll . 
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